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ABSTRACT 

Context. We present multi-wavelength observations of bright quasar HS 1603+3820: the optical data taken with the MMT and Keck 
telescopes, with the 40-50 kms~' resolution, and X-ray data taken by the Chandra X-ray Observatory satellite. 
Aims. The optical spectra contain a very large number of absorption lines from numerous heavy elements. Our goal is to analyse 
those features to obtain constraints on the properties of associated absorbers, to be used in modeling of the quasar intrinsic flux and 
properties of the clouds. 

Methods. We have determined the properties - column densities and redshifts - of the individual components. We derived the X-ray 
properties of HS 1603+3820 and the optical-to-X-ray slope index a ox . 

Results. We found a ox of 1.70, which is on the high end of the typical range for a radio quiet QSO We found 49 individual heavy 
element absorption clouds, which can be grouped into eleven distinct systems. Absorbers from the associated system which likely is 
the one spatially closest to the QSO show large C iv to H i column density ratio, reaching ~20. 

Conclusions. Intrinsic X-ray properties of the quasar are typical. Determination of column densities of ions (including hydrogen) 
gives a strong foundation for modeling of the quasar ionising flux. 

Key words, quasars: absorption lines - quasars: individual (HS 1603+3820) 



1. Introduction 

HS 1603+3820 is a z em = 2.54 quasar discovered during the 
Hamburg/CfA Bright Quasar Survey (Dobrzycki et al. 1996). 
With B = 15.9, it is among the top few brightest quasars known 
in the z em = 2-3 redshift interval (Veron-Cetty & Veron 2006). 
The most striking feature of HS 1603+3820 is, however, the ex- 
treme richness of its heavy element absorption spectrum, and in 
particular the z a bs ~ z em absorption. 

The absorption spectrum of HS 1603+3820 was first de- 
scribed in Dobrzycki et al. (1999, hereafter D99), based 
on medium-resolution (90-100 kms -1 ) spectrum obtained in 
1997 April with the MMT and the Blue Channel spectrograph. 
Although in this spectrum many of the absorption systems were 
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blended, it was clear that the quasar held great promise for the 
analysis of the properties of associated absorbers. 

The absorption systems which are spatially close to the 
quasar can be used as probes for the intrinsic emission of the 
quasar itself, since the conditions in the systems are undoubtedly 
heavily influenced by the QSO flux (e.g. Crenshaw et al. 2003; 
Gabel et al. 2005; Scott et al. 2005, and references therein). 
Systems which contain lines originating from various species 
and various ionisation levels are particularly interesting, since 
they can provide independent way of determining physical prop- 
erties of the absorbing medium, which can cast light at the intrin- 
sic quasar continuum in far UV band as well as the absorbing 
cloud formation. 

With that as the principal goal, we initiated a project of multi- 
wavelength observations of HS 1603+3820, aiming at obtaining 
spectral information at various wavelengths. In this paper we 
present the results of our X-ray and optical observations of the 
quasar. In the forthcoming paper (Rozariska et al. 2007, in prepa- 
ration; Paper II) we use those observations to constrain models 
of the quasar broad band continuum and the cloud confinement. 

In the meantime, Misawa et al. (2003, 2005, hereafter 
referred to as M03 and M05, respectively) observed 
HS 1603+3820 in 2002 March and 2003 July in very high res- 
olution (6.7 kms -1 ), using the Subaru telescope with the High 
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Dispersion Spectrograph. The primary goal of those observa- 
tions was to establish the nature of the absorbers (intervening 
or associated), using the analysis of line time variability, partial 
coverage and profile shapes. High resolution of their observa- 
tions enabled resolving many of the systems into a number of 
narrow components. M03; M05 grouped the absorbers into sev- 
eral systems, designated with the letters from A to H. In our 
optical spectroscopy data we see all those systems and three ad- 
ditional intervening systems. For consistency, and to enable easy 
comparisons, wherever possible we will use those designations 
throughout this paper. 

In a very recent paper, Misawa et al. (2007, hereafter M07) 
expanded the analysis of absorber variability, using new spectra 
from the Subaru and Hobby-Eberly telescopes, and confirmed 
the finding from M05 that C iv in system A is highly variable. 
This paper also described one system (labeled I) which they did 
not have in their older data. This system was seen in our data, 
and for consistency with M07 we rewrote parts of this paper to 
also use this designation. 

Our optical spectra are of lower resolution (FWHM of 40- 
50 kms -1 ), but they complement the Subaru data in two impor- 
tant areas: they have higher signal-to-noise ratio (100-140 per 
pixel in the C iv and Ly a regions) and they cover much larger 
wavelength range (~3700 to ~ 10000 A). The former reduces 
uncertainty in blended regions containing thick lines and helps 
in finding weak lines, while the latter allows analysing lines from 
species inaccessible in the Subaru spectra. The most important 
addition is detection and analysis of hydrogen lines, not seen in 
the Misawa et al.'s spectra, as well as detection of heavy-element 
lines other than C iv in system A. 

This paper is organised as follows. In Section 2 we discuss 
the Chandra X-ray observation of HS 1603+3820. In Section 3 
we present our optical observations of HS 1603+3820, and we 
discuss the individual systems and compare our results with 
those of Misawa et al., discussing the possible sources of dis- 
crepancies. We comment on the absorber variability in Section 4. 
We summarise the results in Section 5. 

2. X-ray observations 

HS 1603+3820 was observed for 9.2 ks with the Chandra X-ray 
Observatory on 2002 November 29 (ObsID 4026), with the 128 
row subarray of the ACIS-S3 chip in VFAINT mode, with frame 
time of 0.4 sec. The effective good time for the observation was 
8.3 ksec. 

The data were reduced using the tools available in the CIAO 
v. 3.3 and SHERPA software packages 1 , using Chandra calibra- 
tion database (CALDB) v. 3.2.1. Because of large calibration 
uncertainties for soft X-rays, the data with E < 0.5 keV were ex- 
cluded from the analysis. We extracted source and background 
spectra from the regions around the quasar using dmextract. In 
total, the net photon yield from the quasar was 145 photons, al- 
lowing for most fundamental spectral fits only. We note that the 
ACIS pile-up is not a problem, both because of the low source 
count rate and because of the low frame time resulting from the 
use of the chip subarray. 

We performed spectral fits with SHERPA, applying the 
newest gain corrections to the event list. We performed the 
fit assuming the intrinsic quasar spectrum to be a power law, 
AE~ y , with fixed Galactic absorption towards HS 1603+3820 of 
N H = 1.3 X 10 21 crrT 2 , as obtained with COLDEN 2 . The best 

1 http://cxc.harvard.edu/ciao/ 

2 http://cxc.harvard.edu/toolkit/colden.jsp 
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Fig. 1. Chandra ACIS-S3 X-ray spectrum of HS 1603+3820 and best fit 
power law plus Galactic absorption. The bottom panels shows the fit 
residuals. 



fit gave the photon spectral index of y = 1.91 + 0.20 and nor- 
malisation at (rest) 1 keV of A = (2.43 + 0.79) x 10~ 4 pho- 
tons cm~ 2 s _1 keV _1 . The slope of the spectrum is within the typ- 
ical range for quasars (e.g. Fiore et al. 1998; Reeves & Turner 
2000). The spectrum and the fit are shown in Figure 1. 

We also attempted fits in which we allowed for additional 
absorption from the quasar, but there was no improvement in the 
fit. 

Combination of the above result with the spectrophotometric 
data from Scott et al. (2000) yields the value of the optical-to- 
X-ray slope, a ox = - log(/ v>2 kev//v,2500A)/ log(v 2ke v/v2500A), of 
1 .70, which is at the high end of the typical parameter range at 
z em ~ 2.5 (Bechtold et al. 2003). 

3. Optical spectroscopy 

HS 1603+3820 has been observed by us with two instruments. 
On 2001 August 19 we used the ESI, Echelle Spectrograph and 
Imager (Sheinis et al. 2002), on the Keck II telescope, with the 
MIT-LL 2048x4096 CCD, with the 0.5" slit. Observing con- 
ditions were excellent throughout the night, with clear sky and 
seeing of the order of 0.7". The combined exposure time was 
4800 sec. The usable part of the spectrum covered the ~3900 
to ~ 10000 A wavelength range and had the resolution of 45- 
50 kms -1 , with the signal-to-noise ratio peaking at ~140 (per 
pixel) near the Ly a emission line, gradually falling to ~ 100 near 
the C iv emission line and then to ~70 in the red end. 

On 2003 June 21 we used the Blue Channel Spectrograph 
on the MMT, in the echelette/cross-dispersed mode, using the 
Lesser 2688x512 "ccd35" detector with the 1" slit. The observ- 
ing conditions were variable. Seeing became worse at the end of 
the night, reaching ~2 arcsec, and the conditions were adversely 
affected by the moon age of ~ 0.75 and smoke from forest fires. 
The combined exposure time was 18,000 sec. The usable spec- 
trum ranged from ~3750 to ~7100 A and had the resolution of 
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40-45 kms~'. The signal-to-noise ratio in this spectrum is in 
the 1 10-130 (per pixel) range throughout the entire wavelength 
range. 

For both sets of observations, the initial reduction was done 
in the standard way with IRAF. We note that despite several at- 
tempts, the extracted MMT spectrum appears to show residual 
flux. Even lines that are clearly saturated and which in the Keck 
spectrum have black cores, do not reach zero intensity in their 
bottom in the MMT spectrum. We elected to not correct man- 
ually for this effect, but we note that this will likely result in 
underestimating the column density of the very strong lines in 
the MMT spectrum. 

After extracting the spectra from individual orders, we per- 
formed the continuum fits to the spectra with the use of the 
FINDSL programme (Aldcroft et al. 1994). The spectra were 
then normalised, rebinned to common wavelength pixels and co- 
added, weighed with the variance. 

Not surprisingly, the trickiest part of this part of the analysis 
were the continuum fits. The spectra are heavily absorbed blue- 
wards of the Ly a emission line. To fit continua in the uncertain 
areas, we applied technique similar to the one described in M03, 
i.e. utilising the spectrum shapes in the unabsorbed orders. We 
note that this effect is less pronounced in our case, since a single 
echelle order in our spectra - which have lower resolution than 
the Subaru spectra - covers a larger wavelength range, providing 
for more possibilities to "anchor" the continuum. 

However, the continuum fits remain uncertain in two partic- 
ularly heavily blended areas. Comparison between continuum 
location in the normalised spectra from Keck and MMT shows 
that there remains a 15-20% uncertainty in those regions. In one 
case it is partly caused by the fact that one of the ESI echelle or- 
ders ends inside the blended region, forcing manual anchoring of 
the continuum fit. In the other case it is the fact that the blended 
region is on the wing of the emission line, which coincides with 
the peak sensitivity of the echelle order; this results in very un- 
certain shape of the continuum exactly where the technique of 
using shapes from unabsorbed regions cannot be applied. We 
will return to this issue in Section 3.1.1. 

We note that our continuum fits generally appear to agree 
well with the Misawa et al.'s analysis of their data. One notable 
exception will be discussed in more detail in Section 3.1.1. 

The next step in the analysis - resolving and identifying ab- 
sorption lines - was a very laborious process, since our spec- 
tra were heavily blended. In order to avoid bias, we performed 
the fits blindly, i.e. we decided to not use the information from 
Misawa et al. As a consequence, our results do differ quite sig- 
nificantly in some cases, especially in the number of individ- 
ual components into which some of the systems were resolved. 
However, since the primary goal of our analysis is to obtain in- 
formation on the the column densities of as many ions as possi- 
ble in order to use them as constraints for modeling of the quasar 
intrinsic emission, the important information is the total column 
density in the system, and not the number of the individual com- 
ponents. In this regard, in the systems we have in common - 
primarily the vicinity of the C iv emission line - our results do 
qualitatively agree with those of M03 and M05. 

The line fits were done using VPGUESS/VPFIT tools. 3 This 
programme fits multiple Voigt profiles, with absorption redshift 
z a bs, column density N and Doppler parameter b as the fitted pa- 
rameters, and it is very well suited for analysis of blended spec- 
tra. However, our spectra are exceptionally crowded and in some 

3 http ://w w w. eso . org/~j li ske/vpgues s/, 
http://www.ast.cam.ac.uk/"'rfc/vpfit.html 



cases we had to use simplifying assumptions in order to prevent 
the programme from performing fits which were obviously non- 
physical. We also had to apply some judgement calls. The cases 
where it was necessary are mentioned in the following subsec- 
tions. We note that the absorption spectrum of HS 1603+3820 is 
so complex and so severely blended that M03 and M05 had to 
resort to similar approach even in the Subaru spectrum, which 
had markedly higher resolution. 

A large part of the work presented in M05 was devoted to 
the analysis of partial coverage of the QSO by the absorbers. It 
is clear that at least some of the absorbers in the vicinity of the 
quasar do not entirely cover the continuum emitting region. In 
our analysis, we accounted for this effect only in the one case 
where it was clearly necessary to obtain physically sound pa- 
rameter fits - the thick C iv lines in system A (see Section 3.1.1), 
which were clearly saturated but which did not reach zero inten- 
sity in the line bottoms. We note that M05 presented arguments 
indicating that other systems - including those with z a b s <k z e m - 
may also show partial coverage. We will discuss this during the 
description of the relevant systems. 

The resolution of our spectra was slightly higher than the ex- 
pected values of the Doppler parameter for heavy element lines. 
High resolution data from Subaru indicated that in many cases 
the metal lines had b as low as ~5 kms -1 , while FWHM of 40- 
50 kms -1 translates (for Gaussian lines) to b a 25 - 30 kms -1 . 
On the other hand, many lines have b comparable to the spec- 
trum resolution or higher. In our analysis of the spectra we thus 
used the iterative approach. We first freed all parameters in all 
fitted lines. Then for all lines for which the resulting value of 
the Doppler parameter was significantly below the resolution, 
we fixed b at the value which corresponded to the value of the 
parameter for other lines for which it was well established (us- 
ing values for other metal lines wherever it was possible and 
hydrogen lines in other cases), scaled with the square root of 
the element's atomic mass, and re-run the fit. In the situation in 
which one or more lines had inconsistent values of the Doppler 
parameter, we investigated the possibility of having the line as a 
blend with other features. 

It is relatively easy to fall into the trap of adding narrow lines 
in order to improve the fits of large features. We avoided that is- 
sue by adding lines only if it resulted in significant improvement 
of the fit, and all cases when this was done were individually 
examined. 

In total, we found 49 individual absorption clouds, which can 
be grouped into eleven distinct "systems." We stress that the sys- 
tems are defined based on their proximity in the spectrum, and 
the fact that two absorbers are considered members of the same 
system does not in principle mean that they they are physically 
related. 

In Tables 1-4 we show line lists from all systems detected 
in the spectra from Keck/ESI. In the tables, the first column 
shows the line designation; the first component refers to the sep- 
arate cloud within the considered system, and the second part, 
after the dot, is the ordinal number of individual line seen in 
this cloud. The second column gives the ion identification (ID). 
The following columns give the absorption redshift, the Doppler 
parameter, column density and (in Table 1 the covering factor, 
Cf, as defined in Hamann et al. (1997). The last column gives 
(where available) the line identification in M05. We mark lines 
for which the derivation of parameters is made uncertain by 
blending and/or by the fact that the line is saturated. The line 
lists from the MMT/Blue Channel spectrum, as well as the en- 
tire Keck and MMT spectra are available on the World Wide 
Web at http://alpha.uwb.edu.pl/mrk/hsl603/. 
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Fig. 2. The fragments of the Keck spectrum of HS 1603+3820 plotted in 
the same redshift range for several absorption features. Systems A-D are 
shown as shaded areas. The Mg n lines are near the very red end of the 
spectrum, and the "emission" features seen in their vicinity are residuals 
of sky lines. The dotted vertical line shows the adopted location of the 
QSO emission redshift. 

In addition, in Table 5 we list the total column densities of 
all ions detected in all systems and compare them to the values 
found in the Subaru data. 

The following section contains more detailed discussion of 
the individual systems. 

3.1. Systems at z a b S ~ z em 

Of the eleven distinct systems, four (A-D) are close to the 
emission redshift of the quasar. The ejection velocities range 
from ~10,000 kms -1 for system A to negative (i.e. z a b s > z em ) 
~1,200 km s _1 for system D. Figure 2 shows the Keck spectrum 
of HS 1603+3820 in the vicinity of the emission redshift of the 
quasar for several absorbers. One can clearly see that some of 
the systems contain lines from several species. 

3.1.1. System A, z abs = 2.44 

Five heavy-element absorbers were recognised in this system in 
the medium-resolution data of D99. At least one of the compo- 
nents, at z a b s = 2.4189 was identified as associated, since the 
Civ(/11548)/Civ(/11550) ratio was 1, indicating that the cloud 



Keck/ESI MMT/Blue Channel 



3000 


i i i i 1 i i i 


i , , , , i , , , 


2000 


1 1 1 1 1 1 1 1 


: i i i i i i i i i 


2500 
2000 
1 1500 

u 

1000 




f 




1500 
1000 




f\ 




500 









500 










5250 5300 5350 5400 5450 5250 5300 5350 5400 

Wavelength [A] Wavelength [A] 



Fig. 3. The sample raw spectra from Keck (left panel) and MMT (right 
panel) as extracted from the echelle order, prior to continuum fitting. 
The dotted lines show the ~5347-5360 A region where our continuum 
fit differs from Misawa et al. See text. 



was optically thick, yet the lines did not reach zero intensity in 
their bottoms, suggesting that the cloud in which the lines origi- 
nated did not obscure the entire continuum emitting source. 

M03 decomposed this system into seven individual compo- 
nents, while M05 increased this number to fifteen systems. Their 
analysis also confirmed that the absorbers do not appear to en- 
tirely cover the quasar, and they found indication that this system 
is variable. These findings were confirmed and expanded in the 
recent paper (M07). 

System A is unambiguously associated with HS 1603+3820, 
and at the same time it has a large relative velocity with respect 
to the quasar (ranging from ~7500 to ~ 10,000 kms~'), which 
means that it is almost certainly located very close to the AGN 
central engine and its properties are very likely determined by 
the radiation from the QSO. This system is thus the ideal candi- 
date for modeling of the quasar continuum. 

In our analysis, we found twenty two components belonging 
to system A. The results are listed in Table 1 ; we show the Ly a 
and C iv lines on Figure 4. For the thick absorbers in the blue end 
of the system we applied the correction for partial covering using 
empirically determined covering factors (listed in Table 1). We 
estimate the uncertainty on the covering factor to be of the order 
of ~10%; the values agree qualitatively with the values found in 
the Subaru data. 

Our results for system A cannot be directly compared with 
the results from Subaru, since at least four of the components 
seen in our spectra (Nos. 19-22 in Table 1) result from the dif- 
ference between the continuum fits in both sets of data. 

In their analysis, Misawa et al. (also in their very recent pa- 
per) treated the wavelength region between ~5345 and ~5360 A 
as unabsorbed (see Fig. 4 in M05 or Fig. 2 in M07). However, 
this region clearly contains absorption in both the Keck and 
MMT data. In Figure 3 we show the raw spectra of this region 
from both Keck and MMT observations. In both cases one can 
clearly see that the continuum does drop bluewards of ~5360 A, 
and we had to fit lines in this region in order to account for 
this drop. This was actually a difficult task, since it was hard to 
resolve the dip into individual systems. We primarily relied on 
the information from the Ly a region, where strong, but heavily 
blended absorption was present. We found that at least four sys- 
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Fig. 4. Absorption lines from System A. Top panel: Civ(/11548), mid- 
dle panel: Crv(/I1551), bottom panel: Lya. Vertical dotted lines with 
numbers show the positions of individual lines as listed in Table 1 . 



tems (i.e. Nos. 19-22) with reasonable properties are necessary 
to satisfactorily fit this region. However, while we are confident 
that the absorbers are present in this region, we have to treat their 
derived physical parameters as somewhat uncertain. We note that 
if the derived densities are taken at face value, then the four sys- 
tems contain of the order 20% of all H i column density in system 
A, but only 2% of the C iv column density; this would suggest 
that the presence of those clouds in the vicinity of System A is 
coincidental and there is no physical relation between them and 
the system absorbers. 

The other place in the system A area where there is a dif- 
ference between our analysis and the results from M05 is the 
~5308-5317 A region. This area contains the Sin(/I1527) ab- 
sorption from system B (it is the feature near A « 5312 A in 
Fig. 4). We obtained a very good fit using eight lines, compared 
with 1 1 lines used in M05. We note, however, that in the earlier 
analysis M03 also used 8 lines. 

Even though it is a heavily blended region, our analysis was 
aided by the fact that our spectra cover large wavelength ranges 
and contain other Si n lines from system B, which we could use 
to constrain the properties of the A 1527 transition. Those lines 
clearly indicated that this line is not as strong as described in 
M05. However, this also indicated that this region had to contain 
additional C iv lines. We added new C iv doublets which are la- 
beled with 5, 6 and 7 in Fig. 4 and Table 1, and we also added line 
No. 3. Both spectra, and especially the Keck spectrum, show that 
those doublets are required. Additionally, we split M05's line no. 
1 into two lines - our doublets Nos. 1 and 2. This significantly 
improved our fit. 

In System A, we also identified the Si iv, Si m, and N v lines. 
We show the Si iv lines in Fig. 5; the lines from the latter two ions 
are in the heavily crowded Ly a forest region. We note that the 
Sirv absorption was not seen in the Misawa et al.'s data, even 
in the spectrum generated by coadding all the Subaru spectra 
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(M07). This is clearly because of high signal-to-noise ratio in 
our spectra. 

We note that the components in system A have a rather large 
value of the C iv/H i column density ratio, higher than 1 in vir- 
tually all cases and reaching ~20 in the extreme case. This is 
seen in both Keck (Table 1) and MMT spectra, and has strong 
implications for the modeling of the intrinsic spectrum of the 
quasar (Paper II). This effect is not seen in other systems, where 
practically all components have Ncn/Nui < 1. 

Our determination of the Civ column density agrees well 
with the value obtained in the Subaru data (Table 5). As men- 
tioned earlier, none of the other ions seen in our spectra was 
seen in the Misawa et al.'s data. 

Finally, we would like to confirm the suspicion first ex- 
pressed in M03 that the feature at 5527 A in the D99 spectrum 
is not a Fen(/I1608) line at z- d bs = 2.4367. This feature is not 
present in our Keck and MMT data, and a close examination of 
the raw 1997 data revealed that it was in fact an artifact caused 
by a cosmic ray hit. 

3.1.2. System B, z abs = 2.48 

Even the medium resolution spectrum from D99 clearly showed 
that this is a strong and complex system. The Ly a absorption 
was very broad, and its non-symmetric profile hinted that it is in 
fact a blend of at least few lines. 

The Subaru spectrum allowed this complex to be resolved 
into 18-19 separate components. As mentioned earlier, our spec- 
tra could be satisfactorily fitted with only eight clouds. We show 
the vicinity of Ly a and C iv lines on Figure 6. The derived prop- 
erties of all lines are listed in Table 2. Note that this system con- 
tains lines from as many as fourteen distinct ions. 

It needs to be noted, however, that in few cases the line 
properties are uncertain. In particular, the seemingly large col- 
umn densities of some of the Si m(/H206) lines (Nos. 2 and 6 
in Table 2) very likely originate from blending with the hydro- 
gen lines. Several other components are also likely affected by 
features in the spectra; some of those features may be unidenti- 
fied lines, while some are very likely to be spectral artifacts. All 
those cases are marked in Table 2. 

Also, direct fit of Sin (7. Sin in Table 2) and 4.Fen in the 
Keck spectrum resulted in unrealistically high column densities. 
That was confirmed by visual inspection. The only workaround 
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Fig. 6. System B: the C iv (upper panel) and Ly a (lower panel) regions. 
Positions of lines from Table 2 are marked. 

was to perform the fit to those lines by trial-and-error, then fix the 
column density for those components, and then refit the whole 
spectrum again. 

The 4.0 1 lines near ~453 1 A are affected by the C iv doublet 
from the J system (see below). 

The reality of the 5.A1 m, 5. Al n and 5.Fe n lines is somewhat 
uncertain. The programme fitted the lines and they were seen in 
visual inspection of the spectrum, and we decided to keep them 
in Table 2. However, we note that the lines themselves and their 
vicinity are rather noisy. 

Total ion column densities seen in our spectra agree well 
with those seen in the Subaru data (Table 5). Differences can 
be seen for some ions, but they can be resolved when blends and 
suspected blends are taken into account. 

3.1.3. System C, z abs = 2.54 

This system (see Figure 7) is very close to the emission red- 
shift of the quasar. In D99 it was considered to be in-falling be- 
cause of the QSO emission redshift of 2.51 adopted in this pa- 
per. However, the arguments presented in M03 favouring higher 
value of the emission redshift, z em = 2.54, are convincing and 
we adopt this value in both this paper and in Paper II. 

The derived properties of C iv absorption (See Table 3) in 
this system agree well with the Subaru data. In our data we also 
see the Ly a line. There is no significant absorption seen in other 
lines, including Mg n. 

3.1.4. System D, z abs = 2.55 

This system is undoubtedly a z abs > z em system (see Figure 8), 
even with the QSO emission redshift of 2.54. The natural expla- 
nation is that this system is infalling onto the quasar. 

Misawa et al. resolved this system into four separate compo- 
nents, and we see the very same absorbers in our data. For all of 
them we identify Ly a and C iv lines. 



Fig. 7. The C iv (upper panel) and Ly a (lower panel) regions of System 
C. Positions of lines detected in the spectra are shown with dotted line 
and labeled with the ID from Table 3. 



The stronger two components (lines No. 3 and 4 in Fig. 8 
and Tab. 3) also show Siiv and Nv lines. We list Nv only un- 
der component #4, but it is very likely a blend of components 
#3 and #4, since both lines in the doublet are broader than what 
would be expected based on the value of the Doppler parame- 
ter of the other lines. The nitrogen absorption is certainly real, 
but the troughs are shallow and are visible only because of high 
S/N of our spectra; resolving it into two components would be 
meaningless. 

In the case of the Keck spectrum, the fits of the Ly a lines #3 
and #4 were diverging when they were combined with lines from 
other ions. We therefore fitted them separately and then fixed the 
values; this is reflected in Table 3, where their values of N and b 
are listed as fixed. 

We also note that the H i line #1 in the Keck spectrum shows 
unusually low value of the Doppler parameter. However, this 
was confirmed by the visual inspection of the spectrum: the line 
is isolated but does appear to be very narrow. Since the corre- 
sponding C iv absorption is unambiguously identified, we have 
no doubts about the reality of this line, and we attribute the un- 
usually low measured value of b to noise fluke. 

3.2. Systems with z abs <s z em 

The spectrum of HS 1603+3820 contains several other heavy- 
element systems at z abs <s z em - We identified seven of them in 
our spectra; the lists of lines are shown in Table 4. Four of those 
systems: E (z abs = 1.89), F (z abs = 1.97), G (z abs = 2.07), and H 
febs = 2.27) were identified in the early Subaru spectra, while 
system I (z abs = 2.18) was described in the recent paper by M07. 

In almost all cases, the properties derived from our spectra 
match reasonably well with those from the Subaru data. The 
two exceptions are the Nos. 3 and 4 C iv lines from system F. 
In the first case M03 and M05 have a doublet of lines of sim- 
ilar strength separated by ~7 kms -1 , which in our spectrum is 
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Fig. 8. As in Figure 7, for System D. 



Fig. 9. The C iv (upper panel) and Si iv (lower panel) regions for System 
E. Markings as in Figure 7. 



not deblended, and the second case is a weak line buried in a 
nearby, stronger line. The total column density derived from our 
data does, however, agree well with the Subaru value (Table 5). 

Although the SirvOU1394,1403) lines in system E (Fig. 9) 
lie in the Ly a forest, they are clearly identified in both Keck 
and MMT spectra. Existence of the second Si iv doublet is rather 
uncertain. We could not see it in the Keck data, although it is 
perhaps present in the MMT spectrum. 

The natural interpretation for the origin of the systems with 
Zabs ^ z em is that they are intervening, i.e. they are caused by ob- 
jects which happen to be placed on the line of sight to the quasar 
but are not related to it, although there are cases in the literature 
where convincing arguments were presented that systems seem- 
ingly far from the quasar could be associated (e.g. Hamann et al. 
1997; Richards et al. 1999). M05 indicated that systems E and 
F could be associated, since their analysis seemed to show that 
some of the components in those systems had partial coverage, 
although the lines did not show any variability. In our analysis 
we obtained satisfactory profile fits without the assumption of 
partial coverage. This accounts for the difference between line 
parameters derived here and in the Subaru spectra. However, the 
total column densities agree well (Table 5). 

In our spectra, we also identify two systems which are not 
seen in the Subaru spectra. We label them J (z a bs - 1 -9266) and 
K (Zabs - 2.2622); the line properties are listed in Table 4. We 
show the two systems in Figure 10. 

System J requires a special comment. The Ly a line for this 
system would be outside of the wavelength region covered by 
both MMT and Keck spectra, and the only identified feature is 
the C iv doublet. Before we proceed, we note that it is somewhat 
surprising that the redder of the two lines, at A Q \, S » 4538 A, is 
apparently not identified in any of the Subaru spectra, since none 
of the line lists seems to contain it. This line is unambiguously 
present in both our spectra. 

What makes this system somewhat uncertain is the fact that 
both those lines coincide very well with the locations of what 
we expect to be the strongest Oi(/I1302) and Sin(/I1304) lines 
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0.0 
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5050 _ 5060 
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Fig. 10. The C iv regions for System J (upper panel) and K (lower 
panel). 



from System B. However, there are indications that C iv is in fact 
present in this location. First, the spectra contain several other 
Si ii lines, and they constrain the column density of that ion quite 
well, suggesting that Sin(/I1304) alone cannot account for all 
of absorption in the redder of the two lines. Second, for Oi to 
be responsible for all of the absorption in the bluer of the two 
lines, at A Q \, S sa 4530 A, it would require Oi column density of 



10 cm , which we consider implausible. 
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The above considerations lead us to believe that system J is 
real, although the fact that both lines are blends makes its prop- 
erties rather uncertain. 

Another of the previously unidentified systems is a system 
at Zabs = 2.2622, which we label K. It is clearly seen in both 
our spectra; obviously, we could identify it because of high S/N 
ratio in our data. The system is displaced by ~400 kms -1 from 
system H. It is thus likely that both those systems originate in 
objects from the same group of galaxies. 

4. Notes on the variability of the absorbers 

Misawa et al. present convincing arguments that the system A 
absorber is variable, but they did not see variability in any of the 
other systems. In general, unlike the Subaru observations, our 
optical data are not well suited for the analysis of absorber vari- 
ability. First, we have data from two epochs only. Second, the 
data were obtained with two different telescopes/instruments, at 
different locations and under different atmospheric conditions - 
not to mention that one of the observations was likely affected 
by smoke from forest fires. Third, one of the extracted spectra 
had issues with residual flux. Fourth, in some places the contin- 
uum fits had to be fudged manually. Fifth, the variability range 
reported in M07 was of the order of Alog/Y ~ 0.3, which is 
only slightly higher than our typical measurement uncertainties. 
All those factors make direct comparison of the derived prop- 
erties of the absorption lines questionable, and especially in the 
regions affected by the continuum fit uncertainty. Any difference 
between absorber properties derived from the two spectra will 
likely contain an unknown component attributable to those is- 
sues. 

Being aware of all those considerations, we examined the 
column densities derived from the MMT and Keck data look- 
ing for any systematic differences between the two datasets. In 
Figure 1 1 we show column densities for the same components 
as derived from the MMT and Keck data. Note that the error 
bars only show formal uncertainties as derived from the fits and 
do not include the unknown uncertainty related to the analysis 
issues. Despite that fact, one can clearly see that the properties 
of components as derived from our two observations agree very 
well. Only very few data points deviate from the A^eck = A^mmt 
line by more than lcr, and the ones that do are practically all in 
the high-TY region, where issues related to the residual flux in the 
MMT spectrum are expected to play a role. 

This conclusion applies to both z a bs ~ z em - including sys- 
tem A, which Misawa et al. have shown to be variable - as well 
as the z a b s «: z em systems. For the latter systems this is not sur- 
prising, although, as mentioned earlier, there were claims in the 
literature (e.g. Hamann et al. 1997; Richards et al. 1999) indicat- 
ing that systems with large redshift difference showed variability 
(see also Hao et al. 2007). 

5. Summary 

The following summarises our findings: 

- The intrinsic X-ray properties of HS 1603+3820 are in the 
typical range. The available X-ray data were insufficient to 
provide meaningful constraints on the properties of the in- 
trinsic absorber. The optical-to-X-ray slope, a ox = 1.70, 
which is at the high end of the typical range at this redshift. 

- In addition to systems seen in D99 and Misawa et al.'s pa- 
pers, we identify two previously unknown intervening sys- 
tems at z a bs = 1.9266, and 2.2622. With eleven complexes 




Fig. 11. Comparison of column densities of all lines for which measure- 
ments exist from both MMT (x-axis) and Keck (y-axis) data. See text 
for discussion. 



of heavy element absorbers, which together contain ~50 (or 
even more) individual clouds - and that not counting Ly a 
forest systems - HS 1603+3820 is truly a spectacular object. 

- Some of the clouds in system A have a high ratio of C iv to 
H i column densities, reaching ~20. Clouds from other sys- 
tems show Ncw/Nm < 1- 

- The data presented in this paper are not well suited for proper 
analysis of line temporal variability. We note, however, that 
the data are consistent with absence of variability for all sys- 
tems. 

- The optical data presented in this paper contain identifi- 
cations and measurements of several lines not seen in the 
Subaru spectra, in particular the hydrogen and magnesium 
lines. 

- Our data contain measurements of column densities of sev- 
eral ions and therefore provide excellent constraints for mod- 
eling of the cloud internal properties and QSO intrinsic flux. 
We will present such modeling in the forthcoming paper. 

Acknowledgements. We thank J. Liske for his help with the VPGUESS/VPFIT 
packages, T. Aldcroft and A. Siemiginowska for helpful comments, and the 
referee for a very careful review. MN would like to thank ESO for hospi- 
tality during his stay there. This work was partially supported by the ESO 
Director General Discretionary Fund, by the National Aeronautics and Space 
Administration through Chandra Award Number 04700527 issued by the 
Chandra X-ray Observatory Center, which is operated by the Smithsonian 
Astrophysical Observatory for and on behalf of the National Aeronautics Space 
Administration under contract NAS8-03060, and by grant 1P03D00829 of the 
Polish State Committee for Scientific Research. 



References 

Aldcroft, T. L., Bechtold, J., & Elvis, M. 1994, ApJS, 93, 1 
Bechtold, J., Siemiginowska, A., Shields, J., et al. 2003, ApJ, 588, 119 
Crenshaw, D. M., Kraemer, S. B., Gabel, J. R., et al. 2003, ApJ, 594, 116 
Dobrzycki, A., Engels, D., & Hagen, H.-J. 1999, A&A, 349, L29 
Dobrzycki, A., Engels, D., Hagen, H.-J., et al. 1996, BAAS, 28, 829 
Fiore, R, Elvis, M., Giommi, P., & Padovani, P. 1998, ApJ, 492, 79 
Gabel, J. R., Kraemer, S. B., Crenshaw, D. M., et al. 2005, ApJ, 631, 741 



A. Dobrzycki et al.: HS 1603+3820 absorption spectrum. I 



9 



Hamann, F., Barlow, T. A., Junkkarinen, V., & Burbidge, E. M. 1997, ApJ, 478, 
80 

Hao, H., Stanek, K. Z., Dobrzycki, A., et al. 2007, ApJ, 659, L99 
Misawa, T., Eracleous, M., Charlton, J. C, & Kashikawa, N. 2007, ApJ, 660, 
152 

Misawa, T., Eracleous, M., Charlton, J. C, & Tajitsu, A. 2005, ApJ, 629, 115 
Misawa, T., Yamada, T., Takada-Hidai, M., et al. 2003, AJ, 125, 1336 
Reeves, J. N. & Turner, M. J. L. 2000, MNRAS, 316, 234 
Richards, G. T., York, D. G., Yanny, B., et al. 1999, ApJ, 513, 576 
Scott, J., Bechtold, J., Dobrzycki, A., & Kulkarni, V. P. 2000, ApJS, 130, 67 
Scott, J. E., Kriss, G. A., Lee, J. C, et al. 2005, ApJ, 634, 193 
Sheinis, A. I., Bolte, M., Epps, H. W., et al. 2002, PASP, 114, 851 
Veron-Cetty, M. P. & Veron, P. 2006, VizieR Online Data Catalog, 7248, 



A. Dobrzycki et al.: HS 1603+3820 absorption spectrum. I 



Table 1. Lines from system A (z abs = 2.44) detected in the Keck spectrum. 



IN o. Line ID 


Zabs 




(km s ) 


log TV 

(cm -2 ) 


Cf 


MUj ID 


l.Hi 


2.41941 ±0.00002 


62.58±6.51 


13.656±0.033 


0.61 




l.CivH 


2.41941 ±0.00002 


72.74±3.42 


14.965±0.044 


0.36 


1;Civ 


l.Nv(| 


2.41941 ±0.00002 


26.74±17.05 


13.235±0.964 






l.Sirv 


2.41941 ±0.00002 


6.50 (fixed) 


12.065±0.182 






2.Hif 


2.42001 ±0.00002 


25.00 (fixed) 


14.496±0.087 


0.61 




2.CivH 


2.42001 ±0.00002 


198.64±21.32 


14.225±0.054 


0.36 




2.Nv|l 


2.42001 ±0.00002 


9.30 (fixed) 


13.008 ±0.257 






2. Si iv 


2.42001 ±0.00002 


6.50 (fixed) 


11.789±0.325 






3.Hi 


2.42344 ±0.00025 


25.00 (fixed) 


13.132±0.037 


0.61 




3.Civ1I 


2.42344 ±0.00025 


61.09±18.07 


13.440±0.068 


0.36 




3.Nv(t 


2.42344 ±0.00025 


9.30 (fixed) 


12.683±0.319 






4.Hi 


2. 42447 ±0.000 19 


25.00 (fixed) 


13.131+0.037 


0.61 




4.Civ 


2.42447±0.00019 


34.72±17.17 


13.090±0.333 


0.36 


2; Civ 


4.Nv jj 


2.42447±0.00019 


9.30 (fixed) 


13.296±0.349 






5.Hi 


2.42554±0.00002 


25.00 (fixed) 


12.420 (fixed) 


0.61 




5. Civ 


2.42554±0.00002 


10.00 (fixed) 


14.313±0.870 


0.36 


3; Civ 


5.Nv(| 


2.42554±0.00002 


9.30 (fixed) 


12.885±0.106 






5. Si m ft 


2.42554±0.00002 


6.50 (fixed) 


12.667±0.041 






6.H1I 


2.42696±0.00003 


25.00 (fixed) 


14.700±0.050 


0.61 




6.Civ 


2.42696±0.00003 


44.13±10.52 


13.322±0.086 


0.36 




6.Nvjt 


(a) 










7.Hi 


2.42833±0.00001 


25.00 (fixed) 


13.941±0.016 






7.Civ 


2.42833±0.00001 


10.00 (fixed) 


13.526±0.055 






7.Nv(t 


2.42833±0.00001 


9.30 (fixed) 


12.650±0.200 






8.H1 


2.43509 ±0.00006 


49.08 ±30.79 


12.361±0.704 






8.Civ 


2.43509 ±0.00006 


38.64±8.37 


13.440±0.058 




4; Civ 


8.Nv 


2.43509±0.00006 


9.30 (fixed) 


12.977±0.477 






8.Sira(( 


2.43509 ±0.00006 


6.50 (fixed) 


<16.000 






9 .Hi 


2.43710±0.00005 


100.56±33.64 


13.091+0.110 






9.Civ 


2.43710±0.00005 


81.40±8.82 


13.970±0.024 




5; Civ 


9.Nv(t 


2.43710±0.00005 


9.30 (fixed) 


13.094±1.220 






9.Sim ft 


2.43710±0.00005 


6.50 (fixed) 


<16.000 






9.Siiv 


2.43710±0.00005 


6.50 (fixed) 


12.050±0.190 






lO.Hi 


2.43753±0.00002 


25.00 (fixed) 


13.262±0.061 






lO.Civ 


2.43753±0.00002 


10.00 (fixed) 


13.596+0.121 




6; Civ 


lO.Nv (( 


2.43753±0.00002 


9.30 (fixed) 


13.717±0.327 






lO.Simjt 


2.43753±0.00002 


6.50 (fixed) 


<16.000 






lO.Sirv 


2.43753±0.00002 


6.50 (fixed) 


11.839±0.292 






ll.Hi 


2.43911 ±0.00030 


138.45 (fixed) 


13.244±0.087 






ll.Crv 


2.43911 ±0.00030 


50.59±40.10 


13.273±0.481 




7; Civ 


ll.Nvjt 


(a) 










lLSimfl 


2.43911 ±0.00030 


6.50 (fixed) 


< 15.376 






12.Hi 


2.43989±0.00019 


25.00 (fixed) 


12.464±0.344 






12.Civ 


2.43989±0.00019 


10.00 (fixed) 


12.825±0.534 




8; Civ 


12.Nvtf 


2.43989±0.00019 


9.30 (fixed) 


12.485±0.880 






12.Sim|) 


2.43989±0.00019 


6.50 (fixed) 


<16.000 






13.Hi 


2.44079±0.00042 


72.42 (fixed) 


13.057±0.103 






13.Civ 


2.44079±0.00042 


77.05±50.41 


13.518+1.886 




9; Civ 


13.Nv(t 


(a) 










13. Si in H 


2.44079±0.00042 


6.50 (fixed) 


<12.734 






14.Hi 


2.44215±0.00030 


25.00 (fixed) 


13.004±0.131 






14.Civ 


2.44215±0.00030 


10.00 (fixed) 


13.787±0.100 




10; Civ 


14.Nv H 


2.44215±0.00030 


9.30 (fixed) 


14.054±0.101 






14.Sim(| 


2.442 15 ±0.00030 


6.50 (fixed) 


<11.587 






15.Hi 


2.44228±0.00021 


8 1.41 ±43.00 


13.106±0.221 






15.Civ 


2.44228 ±0.00021 


123.22±87.76 


13.240±0.305 




11;Civ 


15.Nv(t 


2.44228 ±0.00021 


9.30 (fixed) 


13.808+0.171 






15. Si in (t 


2.44228±0.00021 


6.50 (fixed) 


<11.520 






I6.H1 


2.44325±0.00148 


114.39 (fixed) 


13.086±0.203 
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Table 1. continued. 



Nn I ine ID 




b 


lop N Ct M05 ID 








(cm 2 ) 


16. C iv 


O A A 1 C 1 A AA 1 A O 

2.4432j ±0.00 14s 


1 AA HA SG, r ~A\ 

100. /4 (fixed) 


1 o nnn i a acc n. p „ T 

13./09±0.055 12; Civ 


16. N v tt 


(a) 






16. Si in ft 


2. 44325 ±0.00148 


6.50 (fixed) 


<12.394 


17.Hi 


2.44480±0.00018 


91.30±64.73 


13.093±0.104 


1 / .Civ 


n a a a OA 1 a aaai q 


j\j.\J\)±Z\J. J4 


l j.U/U±U.l /o 14; Civ 


17.Nv ft 


2.44480±0.000l8 


9.30 (fixed) 


13. 094 ±0.433 


17. Si in (t 


2.44480±0.00018 


6.50 (fixed) 


<12.679 


I8.H1 


2.44641±0.00021 


75.27±25.93 


13.154±0.216 


lo.C IV 


O A A £ A 1 1 A AAAO 1 

2.4464 1 ±0.0002 1 


1 AA i CO AA 

109.2 /±5s. 99 


1 O OAT i A I C, P „, 

13.oU/±U.zUz lo; Civ 


18. N v ft 


(a) 






18. Si in (t 


2.4464 1 ±0.00021 


6.50 (fixed) 


<13.483 


19.Hi 


2.44818±0.00002 


25.00 (fixed) 


13. 594 ±0.020 


ly.civ 


O /I A O 1 O i A AAAAO 

2.44s Is ±0.00002 


A 1 C1 i A A iCO 

41. 51 ±49. 68 


1 O AOA i A 

12.920±0.366 


19. N v tt 


(a) 






19. Si in (t 


2.448l8±0.00002 


6.50 (fixed) 


<12.467 


20.Hi 


2.44954±0.00035 


41.27 (fixed) 


12.608 (fixed) 


OA / ' . . ■ 

20. C iv 


O A A A C /I i A AAAO C 

2.44954±0.00035 


62. / /±55. 11 


13.133±0.3ss 


20. Nv ft 


(a) 






20. Si in (t 


2.44954±0.00035 


6.50 (fixed) 


<12.644 


21.Hi 


2.45069±0.00082 


25.00 (fixed) 


12.874 (fixed) 


21. Civ 


2.45069±0.00082 


65.13+108.68 


13.117+0.681 


21.Nvtt 


2.45069±0.00082 


9.30 (fixed) 


11.928+1.265 


21. Si in (t 


2.45069±0.00082 


6.50 (fixed) 


<12.511 


22.Hi1 


2.45174±0.00041 


51.00 (fixed) 


14.245 (fixed) 


22.Civ 


2.45174±0.00041 


53.88±44.25 


13.059±0.340 


22.Nvft 


2.45174±0.00041 


9.30 (fixed) 


12.686±0.923 


22.Siratt 


2.45174±0.00041 


6.50 (fixed) 


<14.411 



(J - heavy-element line in the Ly a forest region of the spectrum; 
\ - saturated line; 

(a) - detected in the MMT spectrum, but not seen in the Keck spectrum. 
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Table 2. Lines from system B (z abs = 2.48) detected in the Keck spectrum. 



No.Line ID 


Zabs 


b 


log A? 


M05 ID 






(kms- 1 ) 


(cm -2 ) 




l.Hif 


2.47723 ±0.00002 


25.00 (fixed) 


14.583±0.228 




l.Crv 


2.47723 ±0.00002 


10.00 (fixed) 


14.223±0.141 


1,2; Civ 


l.Cn 


2.47723 ±0.00002 


10.00 (fixed) 


13.320±0.093 


l.Cn 


l.Nv(| 


2.47723±0.00002 


9.30 (fixed) 


13.617±0.426 




l.Mgn 


2.47723±0.00002 


7.00 (fixed) 


12.123±0.300 




1.A1 in 


2.47723 ±0.00002 


6.70 (fixed) 


11.923±0.469 




l.Aln 


2.47723±0.00002 


6.70 (fixed) 


11.020±1.288 




l.Siiv 


2.47723±0.00002 


6.50 (fixed) 


13.258±0.082 


1,2; Si iv 


l.Simtf 


2.47723±0.00002 


6.50 (fixed) 


13.068±1.434 




l.Sin 


2.47723±0.00002 


6.50 (fixed) 


12.554±0.103 


1,2; Si 11 


2.Hi 


2.47762±0.00004 


25.00 (fixed) 


14.000±0.154 




2.Civ 


2.47762±0.00004 


10.00 (fixed) 


13.643±0.110 


3-5; Civ 


2.Nvtf 


2.47762±0.00004 


9.30 (fixed) 


12.626±0.645 




2.Mg 11 


2.47762±0.00004 


7.00 (fixed) 


11.915+0.416 




2. Aim 


(a) 








2.Ahi 


(a) 








2.Sim(t 


2.47762±0.00004 


6.50 (fixed) 


16.616+0.333 


1; Sim 


2.Sin 


(a) 








3.Hil 


2.47847±0.00001 


25.00 (fixed) 


14.625±0.100 




3.Crv 


2.47847±0.00001 


10.00 (fixed) 


15.265±0.135 


6-10; Civ 


3.Cii 


(a) 






2;Cii 


3.Nv)l 


2.47847±0.00001 


9.30 (fixed) 


13.493±0.072 


1,2; Nv 


3. Aim 


(a) 








3.Aln 


(a) 








3.Shv 


2.47847±0.00001 


26.40±3.23 


12.827±0.077 


3,4; Si iv 


3. Si in (t 


2.47847±0.00001 


26.40±3.23 


14.203±0.206 


2; Sim 


3. Si 11 


(a) 






3,4; Sin 


4.Hi1 


2.47938±0.00001 


25.00 (fixed) 


15.810+1.099 




4.Crv 


2.47938±0.00001 


29.97±1.36 


14.515±0.026 


ll-13;Crv 


4.Cii 


2.47938±0.00001 


29.97±1.36 


14.671 ±0.029 


3,4; C 11 


4.Nv)t,t 


2.47938±0.00001 


38.37±30.95 


13.512±0.087 


3,4,5; Nv 


4.01$ 


2.47938±0.00001 


8.70 (fixed) 


13.754±0.184 


l;Oi 


4.Mgn f 


2.47938±0.00001 


7.00 (fixed) 


15.096±0.159 




4.Mgi 


2.47938±0.00001 


7.00 (fixed) 


11.766±0.209 




4. Aim 


2.47938±0.00001 


6.70 (fixed) 


13.054±0.107 




4.Aln 


2.47938±0.00001 


6.70 (fixed) 


14.518±0.239 


1,2; Aln 


4.Siiv 


2.47938±0.00001 


21.87±5.96 


13.398±0.035 


5,6; Si iv 


4.Simtt 


2.47938±0.00001 


21.87±5.96 


13.810±0.120 


3; Sim 


4.Sin 


2.47938±0.00001 


21.87±5.96 


13.938+0.151 


5,6; Sin 


4.Fera(tt 


(b) 








4.Fenj} 


2.47938±0.00001 


4.65 (fixed) 


13.250 (fixed) 


l;Fen 


5.Hil 


2.47998 ±0.00002 


25.00 (fixed) 


14.110±0.503 




5. Civ 


2.47998 ±0.00002 


10.00 (fixed) 


14.550±0.080 


14,15; Civ 


5.Nvtf 


2.47998 ±0.00002 


9.30 (fixed) 


12.843 ±0.296 


6;Nv 


5. Aim 


2.47998 ±0.00002 


6.70 (fixed) 


11.836±0.573 




5.Aln 


2.47998±0.00002 


6.70 (fixed) 


11.517±0.484 


3; Aln 


5.Siiv 


2.47998±0.00002 


34.32±5.42 


13.075±0.067 


7; Si iv 


5. Si in (J 


2.47998±0.00002 


34.32±5.42 


12.984±0.065 


4; Sim 


5.Sin 


2.47998±0.00002 


34.32±5.42 


11.355+1.345 


7; Sin 


5.Fem (t 


(b) 








5. Fen jt 


2.47998±0.00002 


4.65 (fixed) 


12.843+0.214 




6.H1I 


2.48064 ±0.00002 


25.00 (fixed) 


14.969±0.590 




6.Civ 


2.48064 ±0.00002 


10.00 (fixed) 


14.110±0.076 


16,17; Civ 


6.C11 


2.48064 ±0.00002 


10.00 (fixed) 


13.983±0.053 


6,7; C 11 


6.Nvtf 


2.48064 ±0.00002 


9.30 (fixed) 


12.423±0.371 




6.Mg 11 


2.48064 ±0.00002 


7.00 (fixed) 


12.903±0.115 




6. Aim 


2.48064±0.00002 


6.70 (fixed) 


12.061±0.362 




6.Aln 


2.48064±0.00002 


6.70 (fixed) 


11.968+0.199 


4; Aln 


6.Siiv 


2.48064±0.00002 


6.50 (fixed) 


13.334+0.094 


8,9; Si iv 
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Table 2. continued. 



Nn I ine ID 


•t-abs 


b 


In? N 


M05 ID 









(cm ) 




6. Si in (t,t 


2.48064 ±0.00002 


6.50 (fixed) 


15.883±0.058 


5-7; Si in 


6. Si ii 


2.48064 ±0.00002 


6.50 (fixed) 


13.136+0.122 


8; Sin 


6.Fera ft 


(b) 








7.Hif 


2.48 145 ±0.00002 


25.00 (fixed) 


14.431±0.338 




7. Civ 


2.48 145 ±0.00002 


10.00 (fixed) 


13.979±0.072 


18,19;Crv 


7.Cii 


2.48145±0.00002 


10.00 (fixed) 


13.914±0.052 


8; Cn 


7.Nv ft 


2.48145±0.00002 


9.30 (fixed) 


12.423±0.371 




7.Mg 11 


2.48145±0.00002 


7.00 (fixed) 


13.027±0.112 




7.A1 ill 


^ a 1 A c 1 r\ aaaao 

2.48145±0.00002 


6.70 (fixed) 


12.562±0.150 




7.Aln 


2.48145±0.00002 


6.70 (fixed) 


12.068±0.174 


5; Alii 


7. Si iv 


2.48145±0.00002 


6.50 (fixed) 


13.407 ±0.082 


10; Si iv 


7. Si in ft 


2.48145±0.00002 


6.50 (fixed) 


13.936+0.171 


8; Sim 


7. Si 11 


2.48 145±0. 00002 


6.50 (fixed) 


13.449 (fixed) 


9; Sin 


7.Fem ft 


(b) 








7. Fen ft 


2.48145±0.00002 


4.65 (fixed) 


12.526±0.293 




8.H1I 


2.48226 ±0.00002 


25.00 (fixed) 


15.111+0.100 




8. Civ 


2.48226 ±0.00002 


10.00 (fixed) 


12.537±0.407 




8.C11 1 


2.48226 ±0.00002 


10.00 (fixed) 


14.374±0.072 


9-11; C 11 


8.Nv ft 


(a) 








8.O1 


2.48226 ±0.00001 


8.70 (fixed) 


13 549+0 066 




8.Mgn 


2.48226 ±0.00002 


7.00 (fixed) 


13.014±0.106 




8. Aim 


2.48226±0.00002 


6.70 (fixed) 


11. 751 ±0.670 




8.Aln 


2.48226±0.00002 


6.70 (fixed) 


11.977±0.197 


6,7; Aln 


8.Siiv 


(a) 








8.Simft 


2.48226±0.00002 


6.50 (fixed) 


13.124+0.113 


9,10; Sim 


8. Si 11 


2.48226±0.00002 


6.50 (fixed) 


13.822±0.155 


10,11; Sin 



ft - heavy-element line in the Ly a forest region of the spectrum; 

X - line blended by identified line (0 1 line is blended by C iv from system J); 

t - line is probably blended with an unidentified line; 

\ - saturated line; 

(a) - detected in the MMT spectrum, but not seen in the Keck spectrum; 

(b) - detected in the MMT spectrum, but noise is too strong in order to fit this line in the Keck spectrum. 
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A. Dobrzycki et al.: HS1603+3820 absorption spectrum. I 
Table 3. Lines from systems C and D detected in the Keck spectrum. 



No.Line ID 


Zabs 


b 


\ogN 


M05 ID 






(kms -1 ) 


(cm" 2 ) 




System C z a t, s = 2.54 


l.Hif 


2.53787±0.00003 


25.00 (fixed) 


15.659±0.270 




l.Crv 


2.53787±0.00003 


8.32±4.25 


13.510+0.155 


l;Civ 


System D z abs = 2.55 


l.Hi 


2.55222±0.00001 


25.00 (fixed) 


13.569±0.697 




l.Crv 


2.55222±0.00001 


10.00 (fixed) 


12.585±0.410 


l;Crv 


2.Hi 


2.55383±0.00015 


25.00 (fixed) 


14.414±1.474 




2.Civ 


2.55383±0.00015 


10.00 (fixed) 


13.456±0.790 


2,Civ 


3 .Hif 


2.55409 ±0.00018 


25.00 (fixed) 


16.042±1.916 




3.Civ1 


2.55409 ±0.000 18 


14.45±7.54 


14.407 ±0.863 


3; Civ 


3.Siiv 


2.55409±0.00018 


11. 81 ±8.00 


12.850±0.043 


1; Sirv 


4.Hi<I 


2.55442±0.00005 


25.00 (fixed) 


14.865±1.732 




4.Civ<I 


2.55442±0.00005 


11.40±3.96 


15.200±1.134 


4; Civ 


4.Sirv 


2.55442±0.00005 


10.64±4.40 


13.593±0.085 


2; Si iv 


4.Nv# 


2.55442±0.00005 


38.57±16.44 


13.191+0.105 





f - saturated line; 

# - likely blend of components 3 and 4. 
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Table 4. Lines from systems E-K detected in the Keck spectrum. 
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No. Line ID 


Zabs 


b 

(tans" 1 ) 


logJV 

(cm" 2 ) 


M05 ID 


System Ez abs = 1.89 


l.Hi 
l.Civ 
l.Siiv (t 


(a) 

1.88805+0.00001 
1.88805+0.00001 


10.00 (fixed) 
6.50 (fixed) 


13.873±0.032 
13.036±0.064 


l;Crv 


2.Hi 

2.Civ 

2.Siivtt 


(a) 

1.88851 ±0.00007 

(a) 


10.00 (fixed) 


13.576±0.034 


2; Civ 




System J (new system) z a b s = 


1.93 




1. Hi 

2. Civ$ 


(a) 

1.92658±0.00001 


10.00 (fixed) 


13.727±0.081 




System F z abs = 1 -97 


l.Hi 
l.Civ 


(a) 

1.96399±0.00002 


10.00 (fixed) 


13.699±0.213 


l;Crv 


2.Hi 
2.Civ 


(a) 

1.96420±0.00002 


10.00 (fixed) 


14.185+0.179 


2,3; Civ 


3.Hi 
3. Civ 


(a) 

1.96519±0.00001 


10.00 (fixed) 


14.241+0.030 


4,5; Civ 


4.Hi 
4.Civ 


(a) 

1.96572±0.00001 


10.00 (fixed) 


13.574±0.026 


6; Civ 


System G z abs = 2.07 


l.Hi 
l.Civ 
l.Siiv (t 


(a) 

2.07061 ±0.00005 
2.07061 ±0.00005 


10.00 (fixed) 
6.50 (fixed) 


14.385±0.576 
12.887±0.307 


l;Crv 
1; Siiv 


2.Hi 

2.Civ 

2.Sirv(t 


(a) 

2.07091 ±0.00005 
2.07091 ±0.00005 


10.00 (fixed) 
6.50 (fixed) 


14.487±0.575 
13.491±0.188 


2; Civ 
2; Si iv 


System I z a b s =2.18 


l.Hi 

l.Civ 

l.Siiv 


(c) 

2.17560±0.00002 
2.17560±0.00002 


10.00 (fixed) 
6.50 (fixed) 


13.295±0.036 
12.359±0.069 




2.Hi 

2.Civ 

2.Siiv 


(c) 

2.17606±0.00005 
2.17606±0.00005 


10.00 (fixed) 
6.50 (fixed) 


13.066±0.056 
12.107 ±0.123 




3.Hi 
3. Civ 
3. Si iv 


(c) 

2.17645±0.00001 
2.17645±0.00001 


10.00 (fixed) 
6.50 (fixed) 


13.741±0.033 
12.750±0.040 






System K (new system) z a b s = 


: 2.26 




LHil 
l.Civ 


2.262 14±0.00005 
2.262 14±0.00005 


25.00 (fixed) 
10.00 (fixed) 


14.850±0.615 
12.827±0.332 




System H z a b s = 2.27 


LHil 
l.Civ 
l.Siiv 
1. Si nit! 


2.266 13 ±0.00002 
2.266 13 ±0.00002 
2.266 13 ±0.00002 
2.266 13 ±0.00002 


25.00 (fixed) 
10.00 (fixed) 
6.50 (fixed) 
6.50 (fixed) 


14.075±0.570 
13.200±0.035 
12.527±0.053 
12.784±0.085 


l;Crv 
1; Sim 



D - heavy-element line in Ly a forest region of the spectrum; 

$ - line blended with identified line (C iv in system J is blended with Oi from system B); 
% - saturated line; 

(a) - H 1 lines are not observable in systems E, J, F, I and G, in either MMT or Keck spectra. H 1 lines are seen in system I in the 
MMT spectrum. 
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A. Dobrzycki et al.: HS1603+3820 absorption spectrum. I 
Table 5. Total column densities of ions detected in the Keck spectrum. 



Ion 


log Mot(this paper) log AT t ot(M05) 


logMot(M07) 




(cm -2 ) (cm -2 ) 


(cm- 2 ) 


System A z abs = 2.44 


Hi 


15.124 + 0.030 




Civ 


15.270 + 0.105 15.236 




Nv 


14.527 + 0.102 




Si iv 


12.250 + 0.164 




Sim 


(a) 




System B z abs = 2.48 


H T 
I 1 1 


1 J.7JU It \J. 1 ij 




Civ 


15.471 ± 0.085 15.024 ± 0.047 


15.144 + 0.162 


Cii 


14.956 + 0.025 14.905 + 0.151 


15.570 ± 0.997 


Nv 


14.085 + 0.151 13.663 + 0.275 


13.730 + 0.023 


Oi 


13.965 + 0.116 14.390 + 0.032 


14.608 + 0.028 


Mgii 


15.106 + 0.155 




(0 


13.493 + 0.062 




Mgi 


11.766 + 0.209 




Aim 


13.260 + 0.085 




Aln 


14.522 ± 0.237 12.964 ± 0.023 


13.165 + 0.889 


(if) 


12.542 + 0.113 




Siiv(2a) 


14.037 + 0.032 13.909 + 0.020 


14.017 + 0.014 


Sim 


16.693 + 0.279 14.880 + 0.470 


15.328 + 0.301 


(Hi) 


15.902 + 0.056 




Sin(3) 


14.298 + 0.084 15.977 + 0.471 


15.783 ± 0.354 


Fern 


(b) 




Fen(4) 


13.449 + 0.064 13.600 + 0.030 


13.650 + 0.010 


System C z abs = 2.54 


Hi 


15.659 + 0.270 




Civ 


13.510 + 0.155 13.440 + 0.110 


13.410 ±0.050 


System D z a bs = 2.55 


Hi 


16.081 ± 1.756 




Civ 


15.272 + 0.967 14.919 + 0.200 


16.662 ± 0.060 


Nv 


13.191+0.105 




Sirv(2a) 


13.665 + 0.072 13.514 + 0.091 


13.567 ±0.023 


System E z abs = 1.89 


Civ 


14.050 + 0.024 14.054 + 0.163 


14.034 ± 0.074 


Sirvft 


13.036 + 0.064 




System J (new system) z a b s = 1 .93 


Civ $ 


13.727 + 0.081 




System F z abs = 1 .97 


Civ 


14.618 + 0.072 15.840 + 3.410 


14.505 ± 0.024 


System G z a b s = 2.07 


Civ 


14.740 ± 0.410 14.227 ± 0.054 


14.331 ±0.021 


Sirv(2b) 


13.587 + 0.162 13.260 + 0.054 


13.327 ±0.017 


System I z abs =2.18 


Civ 


13.928 + 0.023 


14.019 ± 0.072 


Si iv 


12.991 +0.032 


12.797 ± 0.015 


Sim 




12.800 ±0.022 


System K (new system) z a b s = 2.26 


Hi 


14.850 + 0.615 




Civ 


12.827 + 0.332 




System H z a b s = 2.27 


Hi 


14.074 ± 0.570 




Civ 


13.200 + 0.035 13.230 ±0.330 


13.914 ±0.141 


Si iv 


12.527 + 0.053 




Sim 


12.784 + 0.085 13.020 + 2.410 


12.807 ±0.119 



M05 data taken from their Tables 1 and 2 (second spectrum in both cases). M07 data taken from their Table 3. 
The symbols denote: 
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(0 - without blended line 4.Mg n in system B; 
(if) - without blended line 4.Aln in system B; 
(Hi) - without blended line 2. Si in in system B; 

(a) - but ion detected, but N measurements are unreliable due to blends; 

(b) - Fe m in system B are detected in the MMT spectrum, but noise is too strong in order to fit this line in the Keck spectrum. 
H - heavy-element line in Ly a forest region of the spectrum; 

%- blend. 

List of Objects 

'HS 1603+3820' on page 1 



